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■ Abstract 

' A quantum evolution model in 2+1 discrete space - time, connected with 3D fundamental 

map R, is investigated. Map R is derived as a map providing a zero curvature of a two 
CO ' dimensional lattice system called "the current system" . In a special case of the local Weyl 

algebra for dynamical variables the map appears to be canonical one and it corresponds to 
known operator- valued R - matrix. The current system is a kind of the linear problem for 
fyr") • 2 + 1 evolution model. A generating function for the integrals of motion for the evolution 

' is derived with a help of the current system. The subject of the paper is rather new, and 

, so the perspectives of further investigations are widely discussed. 
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1 Introduction 

§ ! 1.1 3D integrable models 

In 3D integrable models the Tetrahedron equation (TE) takes place of the Yang - Baxter 
k><( \ equation (YBE) in 2D. Having got a solution of TE, one may hope to construct a 3D 

^ • integrable model. In the case of finite number of states one may construct usual layer - to 

5^ \ - layer transfer matrices T, so that TE provides the commutability of them [Q, ||, ||. Such 

finite states models are interpreted usually as statistical mechanics models. Really only one 
such model still exists, the Zamolodchikov - Bazhanov - Baxter model 0, ||, ^, |^. The 
uniqueness does not mean that 3D world has no interest. 

When 3D R-matrices have infinitely many states, which is more usual in 3D, very natural 
is to investigate a kind of transfer matrices that has no hidden space. We denote such transfer 
matrices as U versus the notation for usual transfer matrix T. Matrices U commute with 
the set of T-s, but have no degrees of freedom when the set of T is fixed. Thus U resemble a 
hamiltonian. Conventionally models with infinitely many states are regarded as field theory 
ones. 

The structure of U is clarified in Fig. |l] for 2D case. Here p and q stand for the spectral 
parameters of the vertices of T, p/q consequently is the argument of the vertices of J7, and 
tjj and a'j are the indices taking values in a finite set. This 2D picture we give just an 
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example for the sake of clearness. The 1 + 1 evolution models, connected with classical or 
quantum bilinear Hirota or Hirota-Miwa equitations on the lattice, are always formulated 
in terms of ?7-type evolution operators, see for example ^ and references therein. 

In 3D, U-matrix appears as an element of a cubic lattice included between two nearest 
inclined planes. We do not draw the graphical representation of 3D U here, we will consider 
sections of the cubic lattice by two, in- and out-, inclined planes mentioned. A two dimen- 
sional lattice appearing in such sections is called the kagome lattice and we will consider 
it in details below. The first who considered U - matrices in 3D, constructed with a help 
of finite - state R - matrix, and constructed some eigenvectors for it, was I. Korepanov 




1.2 3D integrability: usual approach 

So, the origin of 3D integrability is a solution of TE. Those who dealt with it know that 
it is practically impossible to find it directly. For example, even to prove TE analytically 
for an ansatz given and tested numerically is bloody complicated This means, we 

guess, there must be an alternative way of a 3D Boltsmann weights' derivation. 

Primitive way is to find a solution of TE is to consider the intertwining relation for a 
triple sets of 2D L-matrices, 

X! ^^^t^t"* (^ii)i,2 (-^^2)1,3 (■^J3)2,3 = X! i^Z)2,3 (^Z)i,3 (^n)i,2 ^n',h,h ' 

jlj2j3 jl ,32,33 

(1.1) 

where the structure of Li, 2^1, 3-^2, 3 versus ii, 2-^1, 3^2, 3 is the Yang - Baxter structure, and 
the extra indices correspond to the possibility to consider coefficients Ri^f^l^^ as 3D R - 
matrix. TE appears as the admissibility condition for 

-^1,2 -^1,3 ^2,3 il,4 L2A -^3,4 -^3,4 -^2,4 il,4 ^2,3 il,3 -^1,2 ■ (1-2) 

Another scenario is the Zamolodchikov tetrahedral algebra 

*t^2 *^:^3 ^2^3 - E *2:3 ^1^3 K2 ■ (i-^) 

31 ,h ,33 

In the compact form, introducing the formal basis for the indices of R, e(i,j) = \i >< j\, 
RaAc = E ^ut^ ea{juh)eb{j2,k2)e,{n,h) . (1.4) 

3,k 

TE looks like 

Rl,2,3 • Rl,4,5 • ^2,4,6 ' R3,5,6 — ^3,5,6 ' ^2,4,6 ' ^1,4,5 ' ^1,2,3 : (1-5) 

where the alphabetical indices, labelling the numbers of the spaces, are conventionally 
changed to the numerical indices, such change we will make frequently. 

Most amusing thing is that all these really give a 3D i?-matrix: Korepanov's i?-matrix. |^, 
. Korepanov's i?-matrix as well as Hietarinta's one are some special cases of more general, 
complete i?-matrix derived by Sergeev, Mangazeev and Stroganov 1^, and complete R- 
matrix is equivalent to Zamolodchikov - Bazhanov - Baxter's weights in the thermodynamic 
limit. 
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Figure 1: Configuration U • T in 2D. 
1.3 3D integrability: functional approach 

A way to get something else in 3D is to refuse tlie finite number of states in the previous 
approach. Namely, 3D models appear in the local Yang - Baxter equation (LYBE) approach. 
LYBE (i.e. a Yang - Baxter equation with different "spectral" parameters in the left and 
right hand sides) can be adapted to a discrete space - time evolution of the triangulated 
two dimensional oriented surface as a kind of zero curvature condition |l^, |l7| . 

In few words, if a matrix Lij{x), acting as usual in a tensor product of two finite 
dimensional spaces labelled by numbers i and j, with some fixed functional structure and 
depending on a set of parameters x, obeys the equation 



il,2(Xa) il,3(Xb) i2,3(Xc) = -^2,3(x'J il,3(Xb) il,2(xl) , 



(1.6) 



called the local Yang-Baxter equation, so that parameters Xq, x;, and Xc are independent 
and xjj, x'f, and x^ can be restored from (1.6) without any ambiguity. 



X(j /a(Xa,X^,Xc) , X^ /b(Xa,X5,Xc) 

then the functional map R is introduced: 



fc (Xa , Xf, , Xc) , 



R 



a,b.c 



ip{x',,4X) V ^(...) 



(1.7) 



(1.8) 



Due to the difference of the "spectral" parameters in the left and right hand sides of LYBE, 
any shift of a line of a two dimensional lattice, constructed with a help of Li_j(xij), changes 
the set of p ara meters Xij . Partially, if any shift of the lines can be decomposed into primitive 
shifts like (1.6) in different ways, then corresponding different products of R-s coincide. The 
basic example of this is the functional Tetrahedron equation. 
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Suppose we move all the lines of the lattice in some regular way, conserving a structure of 
the lattice. Then the change of parameters be considered as an one-step evolution 

of the dynamical variables x^.j governed by an appropriately defined evolution operator 
U = ritriangles ^- "^^^^ evolution is integrable due to uniqueness of LYBE and The 
partition function for the lattice becomes the natural integral of motion. In terms of the 
transfer matrices, the partition function is the T-type transfer matrix. Being functional, 
these R-operators correspond to something infinitely dimensional. Contrary to the previous 
finite dimensional i?-matrices, there are known a lot of such R-operators. The reader can 
find an interesting set of such simplest functional R-s in ||l9[| . 

A quantization of known functional R-s is still open problem. Simplest functional R-s 
are to be regarded as some functional limits of multivariable R-s with a symplectic structure 
conserving. The problem is to rise known R-s to the complete phase space case, this is done 
just for a couple of R-s. 



1.4 3D integrability: general concept of evolution 

Here we discuss, what else can be invented to get a 3D integrability. 

The main observation is that the relations like tetrahedral Zamolodchikov algebra and 
LYBE have usual graphical interpretation as the equality of the objects assigned to two 
similar graphs. These graphs are the triangles, and we will denote them briefly as A and 
V. Left hand side type graph A corresponds to a product like Li^2 ii.3 -^2,3, and right 
hand side type graph v corresponds to L2,3 L13 Li 2- Algebraic objects are assigned to 
the elements of these graphs. In the case of LYBE these algebraical objects are matrix L 
with the indices assigned to the edges (in the form of subscript, for ii.2 1 and 2 stand for 
the edges), and parameter x assigned to the vertex. From 3D point of view x-s are the 
dynamical variables, whereas L and its indices are auxiliary objects. The equality of l.h.s. 
of LYBE and r.h.s. of LYBE gives the notion of the algebraic equivalence of A and V- 
Note, this form of the algebraic equivalence is not obligatory ! 

3D integrability we can get from any other decent definition of an algebraic 
equivalence. 

In this paper we consider a system, associated with a set of equivalent planar graphs. 
We propose another notion of an algebraic equivalence of equivalent graphs. 

We will deal with all elements of the cw-complex, so we start from recalling the relation- 
ship between the elements of a planar graph and repeating some definitions. 

Consider a graph Gn formed by n straight intersecting lines. The elements of its cw- 

tlItl — 1 ) 

complex are the vertices, the edges and the sites. G„ consists on Ny = vertices, 

Ns = — — ^^^^ — ~^ closed inner sites and Ng — 2n outer open sites, Ne ~ n{n — 2) 

closed inner edges and — 2n outer edges. If two graphs G„ and G'^ have the same 
outer structure, i.e. G'^ can be obtained from G„ by appropriate shift of the lines, then call 
G'„ and G„ equivalent. 

Suppose we assign to the elements of a graph some elementary algebraic (maybe, the 
term "arithmetical" is more exact) objects. These objects are divided into two classes: 
dynamical variables and auxiliary objects (see the interpretation of -/ji,2(x) above in this 
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subsection). Dynamical variables are parameters of graph G„, and auxiliary objects give 
some two dimensional rules of a game (like the summation over all intermediate indices 
in the product of L-s). Dynamical variables plus a rule of game give an algebraic object 
corresponding to the whole graph (like the partition function for L-s). This algebraic (arith- 
metical) object for whole G„ we call the observable object. Denote it 0{Gn)- It depends 
on the set of the dynamical variables. 

Consider now two equivalent graphs, G„ and G'„ . The algebraic problem of the equiva- 
lence arises, 

0(G„) = OiG'J . (1.9) 



If, according to the two dimensional rules of the game, we can get (y_^) choosing the dy- 
namical variables for G'j appropriately for the variables of G„ given, then the algebraic 
equivalence makes a sense. If, m oreover, parameters of G'„ can be restored from the alge- 



braic equivalence condition (L9) without any ambiguity, then this equivalence is decent and 
the integr ability is undoubted. 

The algebraic equivalence usually called zero curvature, and LYBE as the zero curvature 



condition as well as functional evolution models was considered in 16, y_7|. Another 
formulation of the algebraic equivalence, different to the LYBE approach, is Korepanov's 
matrix model (see ^ and references therein). The formulation of the matrix model 
differs from the usual assigning the vertex - type Boltsmann weights to the vertices of a 
lattice, but functional evolution models probably are the same. 
We chose another way. 

The method we use was formulated originally in the classical (i.e. functional) evo- 
lution model was described in |Q , and the quasiclassical case was investigated in . This 
paper contains the overview of the method, and the description of the quantum evolution 
model. The main new result is the generating function for integrals of motion for this 
evolution. 



2 Auxiliary Linear Problem 

In this section we give some rules allowing one to assign an algebraic system to a graph. The 
elements to which we assign something are vertices and sites. First, we give most general 
rules, which do not give an algebraic equivalence of equivalent graphs in general, due to a 
sort of a "gauge ambiguity" . As a special case we find the rules which contain not a gauge 
ambiguity, and so a notion of an algebraic equivalence will be introduced. Then we describe 
the map of the dynamical variables given by the equivalence of 2-simplices, and discuss 
other similar approaches giving this map. 

2.1 General approach 

Choose as a game the following rules: 

• Assign to each oriented vertex V an auxiliary "internal current" (f>. Suppose this 
current produces four "site currents" flowing from the vertex into four adjacent faces, 
and proportional to the internal current with some coefficients a,b,c,d, called the 
dynamical variables, as it is shown in Fig. ^. All this variables, cj) and a,...,d for 
different vertices are independent for a while. 

• Define the complete site current as an algebraic sum of the contributions of vertices 
surrounding this site. 
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Figure 2: The current vertex. 




Figure 3: The Yang-Baxter equivalence. 

• For any closed site of a lattice let its complete current is zero. Such zero relations we 
regard as the linear equations for the internal currents. 

• For any graph Gn the site currents assigned to outer (open) sites we call the "observable 
currents". In part, two equivalent graphs G„ and must have the same observable 
currents - this is the algebraic meaning of the equivalence. 

Clarify these rules on the example of equivalence of G3. As it was mentioned, this is 
usual Yang - Baxter equivalence graphically, A = shown in Fig. |^. Assign to the 
vertices Wj of the left hand side graph A the currents and the dynamical variables 
dtj , bj , Cj , dj , and to the vertices W'^ of the right hand side graph y - the currents (j)'^ and 
the dynamical variables a^ , b^ , c^-, d^. Six currents of outer sites denote as (f>g, and two 

zero valued currents of closed sites - as (j)h and ipa as it is shown in Fig. |^. Then, using the 
rules described above, we obtain the following system of eight linear (with respect to the 
currents) relations: 
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(2.1) 
(2.2) 



(2.3) 



(2.4) 

Given are the currents and the dynamical variables for the left hand side graph. Due to 
0/1 = 0, eq. (2.1), only two currents are independent, let them be 0i and 03. All the 



variables for the right hand side graph we try to restore via the linear system. First, use 
06, 0c and 0d (2.2) to express all 0^. Substitute 0^- into relations for 0e, 0/ and 0g (2.2), 



then it will appear three homogeneous linear relations for two arbitrary 0i and 03, so six 
coefficients of 0i and 03 must vanish. Solving this six equations with respect to the primed 
variables, we obtain 



(2.5) 
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where three polynomials arisen: 



Ai = b3 83 1 ai - ci bi 



A2 = b3 dg ^ C2 d2 ^ 



A, 



32 "^2 d3 C3 



32 d2 



ba C3 



32 b2 ^ dj bj ^ 



Ci di ^ b2 d2 ^ 



C]^ 3j^ 83 C3 



Substituting 0' into 0a = (2^), we obtain the homogeneous linear equation for 



again, and the coefficients of them vanish if 

b'l c'f 1 = K Ai (c2 b^i di br^ + d3 83-1 ai b^^)"^ , 

da 33"' = K A2 (ai b2 + a3 dg ^ C2 d^^)"' , 



(2.6) 



>1, 03 



(2.7) 



C3 b3 



Ka A3 (di ^ 33 C3 ^ 



b2 C2 ^ d3 C3 ^) 



where Aq is arbitrary. The origin of A^ technically is 0a = Aa ■ 4>h- 

This Aa is a sort of a gauge. The origin of it is that due to 0a = we may change it 
0a '—>' ^a'Pa, this givcs Aa > AqAq, Or equivalent 



b'l Aa b'l , d'2 ^ Aa d'2 , C3 Aa C3 



(2.8) 
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Analogous degree of freedom is lost in the map Wi, W2, W3 W{, W'^.W^;. the system of 
the observables is not changed when ^ Xhff'h, i-e. when 



ci ^ Xh ci , 32 1-^ Ah 32 , bs Xh bs , (2.9) 
and the formulae for W': do not change with (2.9). Call such type invariance of the system 



of the observables the site projective invariance (correspondingly, the site ambiguity of 
the dynamical variables) . 

The other obvious invariance (ambiguity) is the vertex projective one. As the conse- 
quence of simple re-scaling of the currents almost nothing changes if 

3i-^3A, bi-^bA, CK^cA, di-^dA (2.10) 

partially in all vertices Wj and with six different Xj and A^ . 

Thus in the most general interpretation: the map Wi,W2, W3 1-^ W[, W2, is defined 
up to projective ambiguity Ai, A2, A3, Xt 1— > A^, A2, A3, Xa- 

Very important feature of all these calculations is that 



we never tried to commute anything 



Return to a general case of graph G„. ANy — 2n{n — 1) free invertible variables 
3i/,by,cy,dy, assigned to the vertices V of G„, we regard as the generators of a body 
B{Gn)- Let Bp{ Gn) be the set of functions on B{Gn) invariant with respect to the vertex 
ambiguity ( 2.10| ). Note in general, for an open graph G„ one may consider B'p(Gn) - set 



of functions invariant with respect to both vertex and closed site ambiguities. But such 
general considerations of Bp for the closed graphs, i.e. the graphs defined on the torus, 
needs a notion of a trace (or of a characteristic polynomials), or equivalent, of an algebra. 
The algebra will be introduced in the subsequent section. 

Consider a little change of G„, so that only one A in G„ transforms into V- Call the 
resulting graph G^. Let the vertices involved into this change are marked as Wi, W2, W3 for 
A and Wl,W2, for y arranged as in Fig. ||. Introduce a functional operator R = Ri,2,3 
making the corresponding map on ;Bp: 

Rl,2,3-'/'(W^l:^2,T^3,...)-R^2,3 = ^W-W^2,W^3:-) , V ^ Bp (2.11) 

where Wj stands for {aj,hj,Cj,dj} forever, and all other vertices except Wi,W2,W3 and 
their variables remain untouched. This R we call the fundamental map. 

Let now G^ be an arbitrary graph equivalent to G„. G^ can be obtained from G„ 
by different sequences of elementary A i— > y in general. Thus the corresponding different 
sequences of R-s must coincide, this is natural admissibility (or associability) condition for 

Gn ^ g;. 

Note that in terms of functional operators the sequence of nai've geometrical transfor- 
mations is antihomomorphic to the sequence of corresponding functional maps. 

The simplest case is the equivalence of two quadrilaterals, G4, and the admissibility 



condition is nothing but the Tetrahedron equation (1.5). And due to the ambiguity of R, 



(2.8,p.S|), any admissibility condition is still equation for Aa-th involved. Note, B'p{Gn) 
introduced previously, is gauge invariant subspace of Bp. R acts on B'p uniquely. Unfor- 
tunately the basis of Bp is not local, and it is simpler to introduce an algebra constraints 
removing the projective ambiguities then to consider B'p formally. 
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Figure 4: Local parameterisation of the vertex 



A way to remove ambiguity from the definition of R, {2^,2/7), is to 



impose some additional conditions for the elements of W, a,b, c, d, such that 



( |2.7| ) would become a consequence of (2^) and the additional conditions 



Complete classification of these additional conditions is still the open problem, and this 
is the main mathematical problem of this approach. 

2.2 Local case: the Weyl algebra 

Here we consider a special local case: suppose first that the elements of two different Wi 
and Wj for the given Gn commute. Destroy also the vertex projective invariance choosing 



aj = 1 for any j forever. Then (2.5) give the expressions for 62,63, d'lc'i ^jdgbg ^, c'i,c'2. 
Suppose also any pair of the variables from W are linearly independent, then 

• the commutability of the elements for different Wj from v gives (after some calcula- 
tions) be = qcb with the same C-number q for any vertex, 

• these relations conserve by the map R, i.e. b'c' = q d b'. 

• Also b^^c^^d appear to be centres, depending on the vertex. 

The gauge ambiguity becomes the ambiguity for these centres. We are going to get a sort 
of quantum theory, b and c are already quantized, so we have to keep all centres to be 
invariant, bj^cj^d^- — b^^^c^^^d^ . This is possible, and further we will threat these centres 
as a kind of spectral parameters. 

Change now notations for the dynamical variables to more conventional, and write down 
the resulting expressions for the map R. New notations for the site currents are shown in 
Fig.|. 

Proposition. • Let the vertex dynamical variables are given by 

a = 1, b = (7^/^ u , c = w, d = KUW, (2.12) 
here u, w obey the local Weyl algebra relation, 

u • w = (J w • u , (2-13) 
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and u and w for different vertices commute, and number k is tlie invariant of tfie vertex, i.e. 
Kjj-, assigned to the intersection of lines i and j, is the same for all equivalent graphs. 

Then the problem of the algebraic equivalence (i.e. equality of the outer currents) of 
two graphs: G with the data 0, u, w, and G' with the data u', w', can be solved without 
any ambiguity with respect to all </>', u', w', and the local Weyl algebra structure for the 
set of u',w' is the consequence of the local Weyl algebra relations for the set of u,w. • 

Write the fundamental simplex map for A = y explicitly. The map R = Ri,2,3 : 



is given by 



where 



R • u, = u^. 



R 



R ■ \Nj = w' • R , j = 1,2,3 





= W2 • A3 , 




= A^i • , 




W2 


= A3 1 • Wi , 




= • U3 , 


(2.15) 


W3 


= A2-1 • ur' , 




= U2 ■ Ai , 





(2.14) 



Ai 


= Ul ^ ■ U3 






A2 


Kl _i 
= — U2 • 




-1 -1 -1/2 ^3 -1 -1 

Ui • W2 q ' U2 • W2 , 

K2 


As 


= ■ W^^ 


- 1'^' U3 


W3 ^ + K3 ■ Us . 



Reverse formulae, giving R ^ , look similar 

Kl 



Ar^ = - 

K2 



K2 



K3 \n[ ^ 



(2.16) 



A2- 



Wo + Ui 



Wo 



q K2 u'2 ■ W2 



(2.17) 



'^3 



• W3 - g 



1/2 ^ 



K2 K2 

The conservation of the Weyl algebra structure 



U3 • W3 + Kl W2 • U3 



Uj ■ Wj = q Wj 



q w,- 



(2.18) 



means that R is the canonical map, hence Ri,2,3 can be regarded as an usual operator de- 
pending on Ui, Wi, U2, W2, U3, W3. The structure of R will be described in the next subsection. 

Now the projective ambiguity is removed, and the current system game gives the unique 
correspondence between the elements of equivalent graphs. This is well defined meaning of 
the algebraic equivalence. Hence all the admissibility conditions (and surely the Tetrahedron 
relation) become trivial consequences of this umambiguity, and we get them gratis ! 
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Mention now a couple of useful limits of our fundamental map Ri.2,3- The first one is 
the limit when ni = K2 = = k, and then k t—^ 0. Denote such limiting procedure via 

Kl = K2 = K3 << 1 . (2-19) 

Corresponding map we denote 3- The conditions for k-s are uniform for whole Tetra- 
hedron relation, 

Kl = 1^2 — = K4 = K5 = Kq « 1 , (2.20) 

so R^' obeys TE. The other case is the limit of Ri,2.3 when 

Ki « K2 = K3 « 1 . (2-21) 

These conditions are uniform for TE again, 

Kl « K2 = K3 << K4 — K5 = Kg << 1 . (2.22) 

Corresponding map we call ri,2,3, and due to the uniformness it also obeys TE. Recall, all 
these maps, R with ki = K2 = K3 = 1, R^' and r were derived previously as a hierarchy 
of R ~ operators solving TE, see |||, |2|, |3|, |2|] . 

2.3 Structure of R 

Remarkable feature of R is its spatial invariance. Change a little the operators on which R 
depends: 

Ti — Kj^^ U2 • Ug"^ ■ ^1 7 r2 ~ ^2 Lll ' *3 ' ■'^S , 1^3 = Kg"^ \N^^ • W2 • A3 . (2.23) 

Then for a, (3, 7 being the cyclic permutations of 1, 2, 3, 

(r/5 • - qVa-Vp) - T^ - ■ {Tfj ■ T^ - qVa- T/j) 

(2.24) 

- (1 -q){i- q^) (r„ - r^) = , 

and 

g r„ • - • r„ - {q^'^ v^-Vp - q-'/^ Vp ■ r,) 

(2.25) 

+ _ ^) (^-1 + qVp - T^) = . 

It resembles SO{3) invariance. 

Give now a realisation of R in terms of more simple functions. First, recall the definition 
and properties of the quantum dilogarithm. Let conventionally 

(x; g)„ = (1 - x) (1 - qx) (1 - q\) ... (1 - q^-'x) . (2.26) 

Then the quantum dilogarithm (by definition) [ p5| , p6| 

^(X) {q''\qU - E i. \ - (2-27) 



and 



00 j^i2 

^(x)- = E 7!:^ ■ (2.28) 
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This function is useful for the rational transformations of the Weyl algebra: 

ij{qx) = (1 - gi/2x)-i yj{x) , yj{q-'x) = (1 - g-^/^x) ^(x) , (2.29) 

hence 

tp{u) • w = w • (1 - g^/^u)"i ■ ■i/;(u) , • u = u • (1 - g^^/^w) • -0(w) . (2.30) 

ip is called the quantum dilogarithm due to the pentagon identity |2^ 

V'(w) • V(u) = tpiu) ■ tPi-q-^/'^ uw) • tPiw) , (2.31) 

this corresponds to Roger's five term relation for the usual dilogarithm. From the other side 
Ip is the quantum exponent due to 

Tpiu) ■ ^p{w) = Tpiu + w) . (2.32) 

Recall, everywhere the Weyl algebra relation u w — (/wu is implied. 

Introduce now a generalised permutation function. Let P(x, y), x • y = y ■ x, is 
defined by the following relations: 



P('?x,y) = y-^P(x,y) = P(x,y)y 
P(x,<7y) = P(x,y)x-i = xP(x,y) 



(2.33) 



and 

For z obeying 
it follows 



P(x,y)2 = 1 . (2.34) 
X ■ z = q^'' z ■ X , y ■ z = q^y z ■ y (2.35) 



P(x,y) ■ z = qf'^ fy z ■ xfy ■ y~^- • P(x,y) . (2.36) 

This function we call the generalised permutation because of usual permutation operator of 
the tensor product is 

P = P(u (g) u"\w® w-i) . (2.37) 
Considering three independent Fq ( 2.23| ), a = 1,2,3, one may see that all them depends 



on three operators U, W and s: 

U = W2"^ ■ W3 , W = wi • Ug"^ - q^l'^ s • U • W"^ = ui ■ u^^ . (2.38) 
U W = (7 W U and s is the center. One can directly verify that 



R = V(k3 U) • V'(W-i) • PiJ— U , ■ W^) . W)-i • 4>{k2 U-1)-i , 



(2.39) 



being substituted into ( 2.14 ), gives (|2. 15 , 2.16 ). On U and W , R acts as follows 



r R . u . R-i = ^ • (w - + u 



K3 



.(W - 9^1/2 ^ s + ACi S • U)"', 

. R ■ W ■ R"i = s • ■ ( W - _^ U ) ■ ( W - q^/^ + ki s ■ U 



-1 



(2.40) 
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When Ki — K2 = = 1, expression (2.39) for R coincides with the operator solution of 



the Tetrahedron equation from ||2^, |2j] . This is the generahsation of the finite dimensional 
3D i?-matrix from = 1 to general q, and the finite dimensional i?-niatrix corresponds to 
the Zamolodchikov-Bazhanov-Baxter model. 

We don't discuss this correspondence here, the reader may find the details concerning 
Zamolodchikov - Bazhanov - Baxter model in |jl|, ^, ^, |l3|, the details concerning the finite 
i?- matrix - in the details concerning the quantum dilogarithm - in original papers 



and operator valued R as the gen eralis ation of finite i? - in l29, 23 , 24 , [30| 



Few words concerning the meaning of ( 2.39 ). All ip-s can be decomposed into the seria 



with respect to their arguments. Substitute these R-s into the Tetrahedron relation (1^) and 
move all the generalized permutations P out. P-s theirself obey the Tetrahedron equation 
and so can be cancelled from TE for R-s. Then twelve ip-s rest in the left hand side of TE, 
and twelve ips rest in the right hand side. The Tetrahedron equation in this case becomes a 
relation resembling the braid group relation in 2D. This twenty-four terms relation can be 
proved directly via the seria decomposition of all 24 quantum dilogarithms. The proof is 
based on several finite q'-re-summations (like g-binomial theorems). This is the first value 



of the formula ( 2.39| ). The second one is that relation ( ^.39 ) gives a nice way to derive the 



finite dimensional complete i?-matrix (just replacing -0-8 and P by their finite dimensional 
counterparts, H, || ^). 

Generalised permutation P(x, y) has no good series realization. Note, if we abolish con- 
dition P^ = 1 for a moment, then formally 

P(x,y) ^ 'Z^^'^x" • (2.41) 

This P obeys P^ = f if one takes the Eulcr definition ^ q" ™ = (5,„,o Hi H. Note 

nez 

that in the manipulations with g-seria the Euler principle "A sum of any infinite series is 
the value of an expression, which expansion gives this series" was never failed. Actually 
P(x, y) is to be defined specially for every realisation of the Weyl algebra. As an example 
mention Kashaev and Faddeev's non invariant realisation of the Weyl generators as shifts 
on the space of appropriately defined functions (p([t]): 

Wj • 'Pm) = [t]j m) , • 0(M) = - [t], ^{[t] : [t], ^ q-' [t],) . (2.42) 

Where [t] is a list of the arguments of if and [t]j is its j-th component. Thus Uj and Wj 
refer to the j-th "pointer" of the list of arguments and hence are not functional operators 
in usual sense. Actually the action of Uj, Wj on V'(M) would be given symbolically by the 
following correspondence: 

ip{[t\) ^ \ti > (g)\t2 > (E>\t3 >«)..., (2.43) 

if the eigenvectors \tj > of the operators Wj might be defined. 
Generalised permutation introduced 

Pi,2,3 = i^l^/^U , s-^W^) = P( /^w^-^wg , -g-i/^ur^wiw^-^swgug-i) (2.44) 

K2 V '*2 
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act of u,, Wj, j = 1, 2, 3, as follows: 



and 



Pi, 2, 3 • wi 

Pi, 2, 3 ■ W2 

Pi, 2, 3 ■ W3 

Pi, 2, 3 • Ui = 

Pi, 2, 3 • U2 = 



' ■ Pi. 



2,3 



/ ^^3 1-. 
W3 • Pi, 2, 3 , 

K2 



(2.45) 



w, • Pi 



2,3 



jj VV2 VV3 



' ■ Pi, 



2,3 



1/2 1-^3 -1 n 

q I . — u^w^ Ug • Pi 

V 1^2 



Pi, 2, 3 • U3 

This gives the following action of P on ip{ti,t2, t^) 
Pi, 2, 3 • v{ti,h,h) = v( 



1/2 -1 n 

Z ' \ — Ui Wi U2 • Pi, 

V '^3 



2,3 



2,3 



(2.46) 



K3 ^3 



t2) 



(2.47) 



where ti, t2, t3 stand on the positions corresponding 1,2,3 of Ri,2,3- 

Another thing to be mentioned is the case of [gj = 1. In this case the quantum dilog- 
arithmic functions should be replaced by Faddeev's integral [ p8| . In few words, it appears 
when one considers the Jacoby imaginary transformation of an argument of ip and q: 



7I/2 



z/6 



-1/2 



Then 



■0f(u) 



and the following expression for '!/'f(u) is valid in the limit of real 6 [g8[ 



V'f(u) (= s{z)) = exp 



sinhTT^ sinliTT^^ ^ 



(2.48) 
(2.49) 

(2.50) 



where the singularity at ^ = is circled from above. 



Return now to map (2.39). The map R conserves four independent operators: 

Wi • W2 , U2 • U3 , s (2-51) 
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and 

H = W]^ • u^^ — g^/^ Ui ■ u^^ • W2 ■ W3 



Ki q Uj^ • W]^ • U2 ^ • Ug ^ + K3 ■ U2 ^ 



K2 q W]^ • W2 ■ U3 ^ ■ Wg ^ + K2 W2 ■ Wg ^ 



(2.52) 

Actually R depends only on two of them, s and H. 
Consider the following product 

a = V^a w^i) ■ u) • ii- c u w) • V(a' w) • u"^) . (2.53) 

Let 

X = a\N~^ +a'w + bu + b'u'^ -q^^^^cuw-q'^/^ab' u~^w~^ . (2.54) 
It is easy to check cr • x = X ' '^^ Hence a as an operator is a function on x- 

a = a{aa , bb' , — \ x) , (2.55) 



a 



I did not find explicit form of function a, only a special case of a when c — b' — 0, then 

V'law"^) V(&u) V(a'w) = ^P{ae-^) tP{a' e) (2.56) 

where 

06*"^ + a' 6 = aw"^ + 6u + a'w. (2-57) 

Nevertheless direct calculations give in terms of a introduced. First, it is convenient to 
rewrite R: 

R = ^(w-i)VX-g^/2AJ3UW-i)P(,/^U,s-iw2)^(-qi/2^U-%)-ii^(^W)-i . 

V 1^2 K3 K3 

(2.58) 

Then 

R^ = N • (2.59) 

where 

N = V(W-^)V'(-q^''^K3UW-^) V'(kiU)V'(s"^W) V(-<7^/^K2s"^U"%) , (2.60) 

and 

D = v(— W)^A(-g'/'^U-iW)V(^U-i)V(— sW-i)V(-<z^/'/cisUW-i) . (2.61) 
K3 K3 K3 K3 

Comparing these with the definition of cr, we obtain 

N = a(s-\K2«3S-\— s|H) , D = a(4s,^s,^s-i|— sH) , (2.62) 

K3 K3 K3 Kl K3 



where H is given by (2.52). 
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X 

Figure 5: Combined vertex 

2.4 Fusion 

One more remarkable feature of the current model is a sort of a fusion. As an example 
consider a planar graph formed by two pairs of the parallel lines. Four vertices arise as the 
intersection points of these two pairs of the lines. This is shown in Fig. ^. 

The vertices are labelled by the pairs of the indices, W\^\, PFi,2, VF2,i and W2,2- Single 
closed site means that there are three independent currents. Let them be the internal 
current 0i.i, assigned to north-west corner, and two currents x and y assigned to southern 
and western semi-strips, x and y are observable currents for this cross considered as an alone 
graph. 

Applying the linear system rules, we obtain step by step 

01.1 = , 

01.2 = y - wi 1 • , (2.63) 

02,2 = W2;2 • a; - 2 ^1.2 ^1,2 ^2^ ' V + '«1,2 ^lA ^1.2 ^1,2 ^2^. ' > 

and from zero value of the closed site current 

02,1 = - W2;2 • X + (Ki 2 Ui,2 Wi 2 w^j - g^/^ Ui 2 w^j) ■ y 

+ (9^^^ Wl.l U1.2 W2;l - 1 1 1 W2"j - Ki^2 Wl.l ^1.2 ^1.2 "^2.1 "^2.2) ' ■ 

(2.64) 

Let further —x' and —y' are the edge variables assigned to the northern and eastern 
semistips. In general they are 

x' = a-x + f3-y + fx-4'^ 

(2.65) 

y' = '^■x + S- y + fy-(j), 
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where 



"2,1 ""2,2 > 



1 
S 



1,2 ""2,1 



,1/2 



-■2,2 ""2,2 



'^1,2 Ui, 2*1^2 ^2.1 ^2,2 : 



1-2,1 "2,1 ""2,2 I 



(2.66) 



2,1 "1,2 "2,1 



1,2 ^1,2 ""1,2 "2,2 ""2,2 
1 



Wo 



'^'1,2 '^'2,1 ^1,2 *1,2 ^2,1^2,2 ' 



and 



h = - q 



1/2 



1,1 "1,2 ^2,1 



"1,1 "1.1 ""1.1 ""2,1 



"1,2 ^1,1 "1,2 ""1,2 ""2,1 
- „l/2 



Ui 



Wo 



f. = 



"2,1 



W 



1,1 "1,2 "2,1 



W. 



2,2 ' 



(2.67) 



"1,1 "-2.1 "1,1 ""1,1 "2,1 



1^1,2 Wi 1 Ui 2 Wi 2 U2,2 W2,2 

-1 



"1,2 '*2,1 ""1,1 "1,2 ""1,2 "2,1 ""2,2 ■ 

Curents x,y,x',y' become the edge currents when we rewrite the cross m Fig. [s] as 
a single vertex with modified (thick) fines; denote it as V7 '--^ {Wi^i, Wi_2, VF2,i, W2,2}- 
Suppose we combine such crosses W (vertices with thick fines) in any way, then zero value 
conditions for the restricted strips (closed thick edges) look very simply: due to the signs 
(— ) in the definition of outgoing x' and y' these zero value conditions becomes "outgoing 
edge current of one thick vertex — incoming edge current of another thick vertex" . Thus 
the strip variables just transfer from one combined (thick) vertex to another, and therefore 
they look like edge variables of tfie thick vertices. 

In the case when for any thick vertex the map x,y i— > x' ,y' (2.65) does not contain extra 
(j), i.e. fx = f y = ill ^-ny sense, then the part of the linear system corresponding to the edge 
variables factorises from the whole current system. If such factorisation exists for a 
graph G then it exists for any equivalent graph G', so sub-manifold of Bp given 
by f^ = f J, = is invariant of a map G ^ G'. 

On this sub- manifold we__can delete all tfie edge currents x = y = ... — 0. In tfiis case 
all corner currents of cross W are proportional to cj) = (j)!.!, and the structure of the "thick" 
vertex becomes the structure of usual vertex. Thus one may define "thick" analogies of u, w 
and K. This phenomenon resembles the usual two-dimensional fusion. 

Write now explicit formulae. Introduce 

1 



Wi 



Wi 



Un 



Wo 



K" 



"l,2"'2,l"'2.2 



k = 



2,1 '*2,2 ""1,1 ""1,2 ""2,1 ""2,2 : 



Wi 1 Wi i Wo 



Wo 



(2.68) 



Ui 1 u 



1,2 "-2,2 "1,1 "1,2 "2,1 "2,2 • 



Uo 1 Uo 



Without mentioning of a representation of the Weyl algebra, its right module etc., suppose 



in (2.65) obeys f^ 



f„ 



0. From this, it follows that 

~-i 



K- 



(2.69) 
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where K is introduced as an "eigenvalue" . On this "subspace" the fusion is defined as 



A(w) = - wi 1 wi 2 , A(u) = 



q 



Jl.l "2,1 I 



A(k) = K 



(2.70) 



The meaning of all these is the following. Consider three "thick" crosses Wj i-^ Wj, 
j = 1,2,3, arranged into "thick" Yang - Baxter - type graphs, A and V- Solving the 
complete problem of the equivalence (with twelve vertices in each hand side) one obtains 
the set of relations like 



3 

E 

fe=i 



Xk ■ fx,. 



Yk ■ fy,k 



(2.71) 



etc., with some Xk and Yk, fx,k and fy^k given by (2.67). A(Wj) we obtain from A(wj) 
applying all eight R-s repeatedly, and A(wj)' is the result of the application of single R in 
terms of A(uj), A(wj) and A{K). 

2.5 Matrix part 



Few words concerning the matrix variables a, /?, 7, S in (2.65). This remark is not important 
for our current approach, but the structure of matrix variables is very interesting. First, 
the map of edge auxiliary variables 



a ■ X + 



(2.72) 



y 



J ■ X + 



f3 -y, 
S ■ y , 

The Yang ~ Baxter equivalence in Ko- 
repanov's interpretation is the Korepanov equation: admissibility of the map of three edge 
variables (x, y, z) assigned to three lines of the Yang - Baxter graph. Let 



appeared in Korepanov's matrix models M, 31 



ai I3i 
71 Si 
1 



X, = 

Then the admissibility is 




^3 = 



1 











as 


/33 





73 


'53 



Xi ■ X2 ■ X^ — x'^ ■ x'2 ■ x[ , 



(2.73) 
(2.74) 



where the primed X-s consist on primed a,/?, 7, (5. Korepanov's equation is equivalent to 
the usual local Yang - Baxter equation for the so-called ferroelectric weights: 



X = 



a (3 
7 5 



L = 



/ 1 \ 
a 
0750 
\ C / 



(2.75) 



where in the numeric case C = a S — P j, and the conventional 2^ x 2^ = 4 x 4 matrix 
form for Yang - Baxter matrix L is used (for the equivalence see |l^ for example). In our 
case the elements of different X-s commute, and the elements of one X obey the algebra 



a ■ P = /3 ■ a , J ■ 6 = S ■ J , 
a ■ "f = q J ■ a , a ■ 5 = q S ■ a 



P-S ^ q S ■ P, 



(2.76) 
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and 



/3-7 = S ■ a — "f ■ f3 . 
From ( 2.76| ) and (2.77) it follows that CP — I^C^Cl — iCj f-iid consequently 

• 7 + 9 7 • /3' - (1 + 9) /3 ■ 7 ■ /3 = 0, 

/3 • 7' + 9 7' • /? - (1 + 9) 7 • • 7 = . 



Hence 



and 



c 



1-g 



(/? • 7 - 7 • /?) , 



(/3 • 7 - g 7 • /?) • 



(2.77) 



(2.78) 



(2.79) 



(2.80) 



Call the algebra of a,P,j,S, given by ( 2.76 ) and ( 2.77 ), as X. Interesting is the following 
Proposition. • Korepanov's equations are nine equation for twelve variables, so Xj are 
defined ambiguously: in general one can't fix one element from a'j^, one fro m d'y , S';^ , and 
one from a'^,5'2. Impose on these three arb itrary elements the simple part of X, ( 2.76 ). Then 
all other r elati ons of X, namely relat ions ( 2.76| ) for the other primed elements and all three 
relations ( 2.77 ) (or, equivalent, ( 2.79 )) for the elements of {X[, Xy, X^}, hold automatically 
as the consequence of Korepanov's equation. • 

This observation, we guess, is a way of a quantization of Korepanov's matrix model. 
Remarkably is that X is the nontrivial algebra. 

As an example consider the case when ^ = 0. Corresponding algebra, Xs=o, contains 
only one nontrivial relation, a • 7 = g 7 • a, and /3 is a center. Being a C - number, f3j 



are to be conserved by the map X, X'. Equations (2.74) contain /32 — Pi P^- Hence /3 



is the pure gauge and one may put /3 = 1. The solution of (2.74) is 





= ("3 + "1 


73)-^ • • "2 


7i 


= / • 7l • "3 


• ("3 + "l • 73)"^ 




= ag + • 


73 




= 7i • 73 




"3 


= «2 • 




/ 

73 


= ("3 + "1 


73) • 71"^ • 72 • • / 



(2.81) 



where / is not fixed by ( 2.74 ), this is the ambiguity mentioned. Permutation relations for the 
combinations of the primed elements, which do not contain /, namely for a'j^, ct^i 72 j "^3 7i 
and 73 7^ , do not contradict the set of the local Weyl algebrae a^- • 7^ — q 7j • . This 
corresponds to the statement of the proposition above. Consider now the Weyl algebrae for 
all primed elements. From this, it follows immediately 



/ 



/ . / • 7j 



Hence / is a C - number, and therefore we may put / 
fixes the ambiguity. 



7, • / . (2.82) 
1 . Thus the conservation of ^5=0 
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The map aj , 



i'j,"fj we've obtained is nothing but ri 2.3, given by the limiting 
procedure (2.21). The identification is a = w^^ and 7 = — q^^^ u • w^^. This case, 

, 1 




X = 



1/2 —1 

— q ' u ■ \N 



(2.83) 



is the quantization of the case (77) from the hst of simple functional maps in p9| . 

The case of general X is rather complicated technically, it is a subject of a separate 
investigation. 



2.6 Co-current system and L-operator 

In this subsection we give another form of the current approach. 

Consider the whole linear system for a graph G defined on a torus (boundary conditions 
assumed). This system is the set of zero equations 

'?^sitc =^ ^vertex ' (/"vertex = , (2.84) 

vertices 

where such equation we write for each site of G, the sum is taken over all vertices surrounded 
this site, and contribution from vertex V, denoted as Wy ■ 0y, is one of ^y, g^^^uy ■ (jjy, 
wycjiv or KyUyWy according to Fig. U and the orientation of V. The toroidal structure 
means that all the sites are closed and the n umber of the sites equals to the number of the 
vertices. Gathering all zero equations ( ^.84| ) together, we obtain the matrix form of them, 

L • $ = , (2.85) 

where we combine the internal currents 4>v into the column <I> and the matrix of the coeffi- 
cients L consists of 

1 , q^^^ Uy , Wy and Ky Uy Wy (2.86) 

for all vertices V of the lattice. L is the square matrix, and explicit form of it depends on 
the geometry G. 



(2.85) can be interpreted as an equation of motion for the action A = $* • L • <&, where 
the row co-current vector $* does not depend on $ and its components are assigned 
to the sites S of G. The equation of motion for $* is • L = 0. Corresponding zero 
equations now are assigned to the vertices, and each such equation connects the site co- 
currents from the sites surrounding this vertex. The problem of the equivalence of A and 
V in terms of co-currents can be formulated as follows: the co-current system for the left 
hand side graph A of Fig. || is 



I'd ■ 



3d • 9 



1/2 



W, 



K2 U2 W2 



<P3 = Vc + ■ (1 ' "3 ■ ""3 ^ Vb ■ 1^3 "3 

and co-current system for the right hand side graph y is 



Ko u, w, = 



= 
= 




(2.87) 



I*/ 

h = 



1/2 ,,/ 



1/2 ,,' 



hho Uq Wq — 



• rt,3 U3 W3 



(2.88) 
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Figure 6: Edge variables. 



The equivalence means that when we remove 0^ from ( 2.87 ) and (/)* from ( 2.88 ), then the 
resulting systems as the systems for 0y are equivalent. 

Consider now co-current equation for single vertex, as in Figs. I or I Let the co-currents 
be 0*, (/)* and where the indices a,b,c,d are arranged as in Fig. ||. The co-current 
equation for this vertex is 







(2.89) 



Suppose we have solved a part of such equations for whole graph G, and obtain (/)*, 0*, (j)^ 
in the form usual for homogeneous linear equations: 



with some multipliers x and y. Then from (2.89) we get 



b*,-q'/^y', or 



6S • X 



where 
with 



U! 



Lu{x, y I u , w ) 



■ y 



y • u 



and y' 



q^/^ U ^ • W + K X ^ 



(2.90) 

(2.91) 

(2.92) 
(2.93) 



Now we may change the interpretation completely. Assign x, y, x', y' to the edges which 
separates corresponding sites. These edge variables are shown in Fig. ||. 

Now we can introduce the auxiliary functional operator L, giving the map x, y i— > x',y', 
as we used to be: 



Lx^y (k, u,w)-x = (jj(x,y|u,w) ^-y - Lx,y (k, u, w) 
Lx^y(K, u, w) • y = x-cj(x,y|u,w) - Lx,y(K, u,w). 



(2.94) 
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With the definition (2.14), L operators obey 



Ly,z(K3, U3, W3) • Lx,z(k2, U2, W2) ■ Lx,y(Ki, ui, wi) • Ri,2,3 = 

(2.95) 

= Rl,2,3 ■ Lx,y(Kl, Ul, Wl) • Lx,z('«2U2, W2) • Ly,z(K3, U3, W3) 

Moreover, Local Yang-Baxter relation 

Ly,z(K3, U3, W3) • Lx,z('«2, U2, W2) • Lx,y(Kl, Ui, Wi) = 

(2.96) 

= Lx,y(Kl,Ui,Wi) • Lx,z(k2,U2,W2) ■ Ly^z (k3 , U3 , W3) 

as a set of relations for uj,,wj., with Ufe,Wfe given and with x, y, z arbitrary, gives again the 
map ( ^.15| ) uniquely ! Thus the kind of the local Yang - Baxter relation appears and for 
our current approach. 

Conclude this section by few remarks concerning the functional maps. All the maps 
introduced are connected to several graphical manipulations. Usually we combine such 
manipulations (A i— > \y oi x,y i-^ x',y' etc.), and write the sequence of the dynamical 
variables' sets obtained S i-^ E', in the direct form 

E = So ^ Si E2 • • ■ S„_i E„ , (2.97) 

where Aj stands for j-th manipulation, which allows us to calculate in terms of previous 
variables Ej_i. The same result. Eg ^ E„, can be obtained as 

Ai A2 ... A„ • Eo = E„ • Ai A2 ... An , (2.98) 

where A^- is a functional operator corresponding the manipulation Aj. Remarkable is the 
reverse order of the operators with respect to the nai've manipulations. Note that the direct 
order we obtain considering the "pointer" action of the operators, as it was mentioned in 
the previous subsection, but the "pointer" action is not suitable for the quantization. 



3 Evolution system 

In this section we apply operator R defined in the previous section to construct an evolution 
model explicitly. Due to the current system's background we formulate this model in terms 
of the regular lattice defined on the torus, its motion, its current system and so on. 

The main result of our paper is the generating function for the integrals of motion for 
the evolution. The derivation of the integrals is based on the auxiliary linear problem. 

3.1 Kagome lattice on the torus 

An example of a regular lattice which contains both A and y - type triangles is so-called 
kagome lattice. As it was mentioned in the introduction, the kagome lattices appear in the 
sections of the regular 3D cubic lattices by inclined planes. Thus the kagome lattice and its 
evolution corresponds actually to the rectangular 3D lattice and thus is quite natural. The 
kagome lattice consists on three sets of parallel lines, usual situation shown in Fig. |^. The 
sites of the lattice are both A and v triangles, and hexagons. 

For given lattice introduce the labelling for the vertices. Mark the A triangles by the 
point notation P, and let a and b are the multiplicative shifts in the northern and eastern 
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Figure 7: The kagome lattice. 



directions, so that the elementary shift in the south-east direction is c = a ^b. Nearest to 
triangle P are triangles aP, bP, cP, a^^P, b^^P and c^^P. Some of them are shown in Fig 

0- . . 

For three vertices surrounding the A-type triangle P introduce the notations (1,P), 

(2,P) and (3,P). These notations we will use as the subscripts for everything assigned to 
the vertices. 

This kagome lattice we define on the torus of size M, formally this means the following 
equivalence: 

P ^ h"' P ^ c^^ P ^ P . (3.1) 

Since the notion of the equivalence, we may consider the shifts af all inclined lines through 
the rectangular vertices into north-eastern direction as it is shown in Fig. ^. It is easy to 
see that Fig. ^ is equivalent to Fig. |[ The structure of the kagome lattice conserves by 
such shifts being made simultaneously for all A-s, but the marking of the vertices changes 
a little. This is visible in Fig. ||. 

Give now pure algebraic definition of the evolution. The phase space of the system is 
the set of 3 Weyl pairs Uj.p and Wj_p, j = 1, 2, 3, P = a" 6'^ Pq, where Pq is some frame 
of the reference's distinguished point, and the toroidal boundary conditions mean 

(3.2) 

Wj.aAf p = WjM'-' P = ^j,P ■ 



The phase space is quantized by the definition. Let Uj p, p for P fixed are given by ( 2.15| ), 
so that the map {uj p, Wj,p} ^ {u^ p, p} is given by the operator 

7^ = n (3-3) 
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Figure 8: Geometrical representation of evolution. 



where Rpi acts trivially on the variables of any triangle P ^ P' . Note, we suppose Kj^p do 
not depend on F, 

Kj- p = Kj , (3.4) 

so that with respect to k-s the translation invariance of the lattice is assumed. Define the 
superscript as follows: 







W*p 


* 

2,aP 






3,6P 


= U^j p , 


W3,6P 



n,p ' 



"2,P I 



"3.P 



(3.5) 



This identification means following: u*p,w*p are the variables which appear on the places 
of previous Uj,p,Wj_p according to Fig. ||. The evolution operator U : {uj_p,Wj_p} ^ 
{u*p,w*p} we define as usual: 



U • p ■ U"^ = u*p 



U • p • U ^ = w 



(3.6) 



Regard the primary variables {uj^p,Wj^p} of the given lattice as the initial data for the 
discrete time evolution, 



Uj.p = Uj^p(O) , Wj,p 
The evolution from < = n to < = n + 1 is just 



Wj,p(0) 



Uj^p(n+1) = U • Uj,p{n) ■ U 



Vj,p(n+ 1) 



U • Wo p(n) • U 



(3.7) 
(3.8) 



Surely, the map U is the canonical map for the Weyl algebrae, so that U is the quantum 
evolution operator. Further we'll consider mainly the situation for t = and the map 
from t = to t — 1. We will omit the time variable and write / instead of /(O) and 
/* = U • / ■ U^^ instead of /(I) for any object /. Due to the homogeneity of evolution 
(3.8, p.6| ,3.5) our considerations appear to be valid for a situation with t — n and the map 
from t = n to t==n+l. 
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3.2 Linear system 

Investigate now the linear system for the quantum system obtained. 

Assign to the vertex (j, P) of the primary (t = 0) kagome lattice the internal current 
(pj^p. The linear system is the set of 3M^ linear homogeneous equation for SAf^ internal 
currents 



r fi,p 

h,p 
h,p 



-,1/2, 



1 /2 

Wl,p • (/)l,p + cj)2,P + q ' U3^p • (/)3_p = 

g^'-Ui,p ■ (?il,P + K2U2,aPW2,aP ' 02,aP + Ws^fep • 03, ftp = 

01,a-ip + '«lUl,6-ipWi b-l p • ^i^b-ip + W2,p ■ 02, P 
1 /2 

+ 9 Ufc,b-ip • 02,b-ip + 03,a-ip + K3U3,pW3_p ■ 03_p = 



(3.9) 



Here we have introduced absolutely unessential notations /j.p just in order to distinguish 
these equations, fj^p are assigned to the sites. Due to the homogeneity we may impose the 
quasiperiodical boundary conditions for 4>i.p: 



4>j.a^' P — A (t)j,p , 0j,hM p 



B 



(3.10) 



It is useful to rewrite this system in the matrix form as (2.85), F = L • $ = 0. 
First combine 0j,p with the same j into the column vector with components, so as 
($^)p = 4'j,P- Introduce matrices Ta and T;, as 



{Ta ■ $,)P 



•,aP , (Tb ■ ^j)p — 4'j,bP 



Due to (|3.10|) 



t5 = A , T'" = B 



M 



(3.11) 



(3.12) 



Combine further Uj.p and Wj^p with the same j into diagonal matrices Uj and Wj with the 
same ordering of P as in the definition of , 



Uj ~ diagp Uj^p , Wj = diagp Wj,p 



Obviously, 



(Ta ■ ■ T-i)p = u,-,p , (T, . u, ■ T-')p = u,,6P , 
and the same for w,-. 



(3.13) 
(3.14) 



Combine further <I>i, <I>i, $3 into 2>M^ column Then from (3^) the matrix L can be 
extracted in the 3x3 x A-P block form: 



/ wi , 1 , U3 \ 

g^/^ Ui , Ta K2 U2 W2 , Tf, Wg 

V + Ui Wi , W2 + U2 , + K3 U3 W3 / 



(3.15) 



Recall, system L • $ ~ is 3M^ equations for 3M^ components of $. 

Introduce now co-currents. As it was mentioned, L • $ = we regard as the equations 
of motion for 2D system with the action 



^ = $* . L • $ . 



(3.16) 
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Figure 9: Co-currents on the lattice 

The block form of the co-currents $* is thus fixed from the form of L, or from (^jj). 
Equations of motion for are F* = $* • L = 0, and in the component form 

fl.P = 01,P • Q^^^ p + (l)2.b-^P + 02,a-ip • «1 Ul,p Wi p + (I>I p ■ Wi p , 
< /2*P = <l>tp ■ 1^2 U2,p W2^p + 02,P • Q^^^ ^2,P + 4'*2fi-^P ' W2,p + 03,aP > (3-17) 
. /a.P = </'l.P • W3,p + 02.P + '/'2,a-ip • U3 p W3 p + (ji^j^p ■ q^l"^ U^ p . 

Here ]* p corresponds to (j, P)-th vertex. The assignment of the co-currents is shown in 
Fig. I ^ 

Elements of F* — $* • L have the following remarkable feature: coefficients in /*p 
belong to the algebra of u^ p, w^ p only. We will use this in the next subsection. 

3.3 Properties of L and the quantum determinant 

Consider first the general properties of equation $* • V = for a matrix V similar to L 
introduced: 

V = ||v,,fc|| , (3.18) 

with the commutative columns, 

Vj,/ : Vj,fc • ^i',k' - ^j'M' ■ '^j.k = if k' ^k. (3.19) 

Such matrices have the following properties. 
Property 1: Consider a system 

Zj ■ Vj.k = Oik (3.20) 

3 

with Qffc being C - numbers, as the system for Zj. Then for k ^ k' 

ak ak' - ak' "fc = ^ Zj> ak yj',k' - ^ zj Uk- \/j,k = ^ {zj'Zj-ZjZji) Vj^k '^j',k' = . 
j' j id' 

(3.21) 
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Matrix ||vj j. • v^/ is non-degenerative in general, so the last equality gives immediately 

■ ^j' ~ '^■y ■ ■ (3.22) 
Consequence: Let ||vi j|| is the inverse to ||vj.fe|| matrix: 

• Vj^fc = ^ Vij • Vj^fc = , (3.23) 

d j 

then 

V? • Vijv - Vjjv • Vij- = . (3.24) 

Property 2: Because of in ||vj^fc|| non-commutative elements belong to the same column, 
the algebraic supplements Vk,i as well as the quantum determinant det (v) are well defined. 
Here we've used the notation "det" as the formal operator-valued determinant 

det I |v,,, 1 1 = ^(-1)'^ n ^^-.-0) ■ (3.25) 

Vi.j and det (v) are polynomials of Vj^fe such that in each summand all multipliers belong to 
different columns and thus commute. Moreover, if in ||v|| two rows coincide, then det (v) = 
0. Hence 

yj,k ■ Vk,i = (5,,; det (v) . (3.26) 

k 

Note, yj.k ■ Vk,i = Vk,i ■ \)j,k- 

As it was mentioned previously, sometimes it is useful to introduce formally a module 
for the body of ||vj_fc||. Here we do this, introducing (p* and (j>Q which belong to such formal 
module. This allows us to formulate the following 
Consequence: Consider now the system of co- vector equations 

($* . V), = ^ 0* • V,, = . (3.27) 

j 

Due to property 2 all (p* belong to the null space of det (v): 

(/}* ■ det (v) = . (3.28) 

From the other hand side, (/)*-s are connected by Zjj' - some rational functions of Vj^k'- 

cP*. = 0* • z,,,. . (3.29) 

Property 1 provides the commutability of Zjjv, hence a solution of ( ^.27 ) can be written as 

(j)* — (j>Q ■ Zj , Zj ■ Zji — Zji ■ Zj , 00 ■ "^st (v) = , Zj ■ det (v) = det (v) • Zj , (3.30) 

where in general Zj ^ Zj. 

Apply now both properties and their consequences to L given by ( 3.15|) . First, for any 
representation of the Weyl algebrae the null subspace (j)Q of whole Gilbert space is defined, 

(/)o ■ det (L) = . (3.31) 

The existence of 0q means the solvability of $* • L — 0. Corresponding Zj have the lattice 
structure, Zj^p. These commutative elements are assigned to the sites of the kagome lattice, 
and observable are Zj p z~,^p,. These operators connect the co-currents in different sites, 
and thus Zj^p actually give the realisation of the path group on the kagome lattice. 

Another important thing is that due to T*^ — A and T^'^ = B, det(L) is a Laurent 
polynomial with respect to the quasimomcnta A and B. 
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3.4 Evolution of the co-currents and integrals of motion 



Consider now the shift of the inclined Hnes giving the evolution. The internal currents as 
well as the co-currents change, and we can trace these changes. 
Introduce two extra matrices, K and M: 



where 



Ao 





( " 


, Ao 


K = 










1 1 




— q 




1 w^^i 


q 


-1/2 


A, + ^ 

K2 



\ 



J 



, -1 -1 

+ — Ui W2 U3 , 

Ai + T-i ^ A3 

K2 



(3.32) 

(3.33) 
(3.34) 



with Aj standing for the diagonal matrices with the entries given by( 2.15| ) correspondingly, 
and 



/ 



, U'2 Ta 



M 



'^1 -1 -1 / / 

— W2 U2 Wj , 



K3 
K2 



, — W2 1 U2 1 U3 W3 Tft 



V 



-1 -\ I rr 

, W3 W2 Ta , 



(3.35) 



/ 



Apply the evolution operator U to L: L* = U • L ■ U ^ . 
( 



"1 > 



1 



^1/2^, 



'«2U2W2T'a , 



(3.36) 

The following relation can be verified directly: 



K 



L M 



(3.37) 



M in general is the matrix making p ^ 4'k,P: and K makes (f)^ p 4>^p- Also K and 
M admit 

Kh^K + L- N, Mh^M + N-L* (3.38) 

with arbitrary N. One can prove the following 
Proposition: • K • det(L) = det(L) • K • 

One can understand this in other terms. Since $* • L = can be solved for t = 0, then 
for t = 1 equation $** • = must also be solved because they are bounded by simple 
linear relations. Hence subspace (/)o must coincide with ('/>o)*, i.e. 



det(L*) = det(L) • D 



(3.39) 
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with some operator D. One may hope, D is not too comphcated, and ( 3.39| ) is not trivial. 

Careful analysis of K and M shows that this D does not depend on the quasimomenta 
A and B. In the functional limit q^^^ i— > ±1 one may easily calculate the determinants of K 
and M, both them are proportional to B^^ , and this term cancels from the determinants 
of the left and right hand sides of (3.37). This is so and in the quantum case. 

Hence D in (3.39) is a ratio of any A,B- monomials from det(L) and det(L*). Element 
D can be extracted, say, from A^'^ B^^' component of det(L): 

^ = n ^r.p ■ n <p ■ (3.40) 

p p 

This means that we can introduce a simple operator d: 

D ^ d ■ d*-^ . (3.41) 

Thus 



is the invariant of the evolution, J* = J, i.e 



J = det ( L ) • d (3.42) 
= J, i.e. 

U • J = J • U . (3.43) 

Decompose J as a series of A and B, 

J = ^'^ -""./^ ' (3.44) 

a,/3en 

where a and j3 are integers and their domain (Newton's polygon) 11 is defined by |q;| < M, 
|/3| < M and |a + /3| < M . Quasimomenta A and B are arbitrary C - numbers, and the 
invariance of J means the invariance of each Sa^p- From the other side, J is a functional of 
the dynamical variables of the lattice, i.e. 

Ja,/3 = Ja,/3({Uj,p , Wj^_p}) . (3.45) 

Surely, due to the homogeneity of the lattice these functionals does not depend on time 
layer, and hence the conservation of J, J* = J, means 

Ja,/3({u,,P,W,-p}) = J„,^({u*p,W*p}) , (3.46) 

i.e. functionals ^a,fj give the integrals of motion in usual sense. 
Note further, 

$* . K ~ $* • (3.47) 

where it is suppose d 0^ • U^^ ~ 0q, and ( |3.47 ) gives the linear action of U on Zj^p. To get 
the equality from ( ^.47 ), one has to normalize only one component of $*. 

Some elements of det(L), corresponding to the border of the Newton polygon 11 of 
J(A, B), can be easily calculated. Operator d introduced is defined up to any integral of 
motion. The simplest integrals are 

jl = W U2,P U3,p , j2 = ]^ Ui^p W3^p , is = n *2,P , (3.48) 

P P P 
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and the convenient choice of d is 

d 



d can be absorbed into det, 
where 



= det (L(°)) , 



(3.49) 



(3.50) 



,(") - 



7-1/2 ^-1 



7I/2 



K2 W2 



q 



T/, u 



b U3 



\ + Tf,^ Hi Wi 



71/2 



U2 ^ W2 + Tf, 1 , 1 W3 1 U3 1 + K3 / 



(3.51) 

Whole number of Ja,0 is 3Af ^ + 3Af + 1, and there are + 1 independent between them, 
and between these one can choose only 3M2 commutative, so J gives the complete set of 
integrals. (As to whole number of summands in J, e. g. for M = 2 it is 1536 = 2^ 3.) The 
existence of 3 abelian integrals is the hypothesis tested for small Ms. 

All integrals corresponding to the boundary of domain n, \a\ = M, |/3| = M, \a+P\ = M, 
are equivalent to the following 3M elements; 



n < 

a 

n ^ 



2,ai + ''b-'' Pa "^.^ai+^b-" Po ' 



(3.52) 



n ^1 



aob'Po '^a.a^b'^Po ' 



where Pq is some frame of reference's point as previously. Note, Vj are not Ta,Tb ~ invariant, 
but restoring this invariance in any way, one obtains the invariants of U. Between Wj,TJj, Vj 
one may choose 3M — 1 commutative elements. Inner part of 11 gives — 3M + 1 highly 
complicated independent integrals, which gives g = 3Af ^ — 3Af + 1 commutative (up to 
(3.52)) independent elements. Note, g is the formal genus of the curve J(v4, B) = const. 



3.5 Walks on the lattice and the integrals of motion 

Give now a geometrical interpretation of the integrals of motion. This interpretation follows 
directly from the analysis of the determinant. Every integral of motion is a sum of monomials 
associated with walks on the lattice such that all the walks have the same homotopy class 
with respect to the torus on which the kagomc lattice is defined. 

It is useful to formulate the walks in terms of general vertex variables a, b, c and d as 
in Fig. ^. Recall the shorter notation W = {a,b,c, d} for the dynamical variables' set. 
Consider the matrix L in this general case. Each row in L corresponds to a vertex of the 
lattice, and each column of L corresponds to a polygon (i.e. to a site) of the lattice. Thus 
det(L) consists on the monomials, each of them corresponds (up to a sign) to a product of 
different Wj^p such that: 
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Figure 10: Fixed outlets for the lattice walks 



• for any vertex {j,P) only one of aj^p, ^j,p, Cj,p, dj,p is taken in this monomial, and 

• for any site only one of surrounding a, d is taken in this monomial. 

Take the lattice and mark the places of the vertex variables a, d, corresponding to the 
monomial, by the arrows, ingoing to the corresponding vertices. Thus for any site and for 
any vertex we have painted only one arrow. 

In order to get a purely invariant functional, we have to multiply det(L) by an integrating 

monomial, in general case this monomial is, for example, Y\ ^Tp "^2 p '^3 p- This choice of 

p 

the integrating multiplier corresponds to element d given by ( ^l9| ). It is easy to see that 
this monomial has the same structure as described above. But due to the power —1 we may 
interpret geometrically this monomial as the set of outgoing arrows. 

The system of the outgoing arrows is thus fixed and shown in Fig. |l^ for each A - type 
triangle of the lattice. For the system of the outgoing arrows and any system of ingoing 
arrows the following is valid: 

• for any site there exist exactly one outgoing arrow and exactly one ingoing arrow, and 
they may touch the same vertex, and 

• for any vertex there exist exactly one outgoing arrow and exactly one ingoing arrow, 
and they may belong to the same site. 

Hence there is the unique way to connect all the arrows inside each site so that a walk 
appears. 

So, the walks we consider, obey the following demands: 

• the system of outlets of the walk is fixed and given by Fig. |l^, 

• the walk visits any site only once, 

• the walk must visit all the sites and 

• the walk must visit all the vertices. 

For any walk W let a(yV) be the number of the components of the connectedness (i.e. the 
number simply connected subwalks). 
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Let now walk W belongs to a given homotopy class a A + (3 B of the torus, where cycle 
A corresponds to T^^ and cycle B corresponds to T*^, and denote such walk as yVa,f3- 

To a walk given assign a monomial according to the following rules: let the walk pass 
through vertex {j, P) so that the walk ingoes the vertex form the side x S W^j,p, and outgoes 
the vertex from the side y G Wj^p. Then the multiplier corresponding to {j,P) is x • y~^. 
The monomial Jw is the product of such multipliers corresponding to all the vertices. Thus 
the reader may see that each monomial we construct gains the structure of an element of Bp, 
described in section "Auxiliary linear problem" , subsection "General approach" : monomial 
Jw, 

Jw = •■• X ■ • x' • y'-i •■• , (3.53) 

X and y are assigned to a same vertex, so x • y^^ does not contain the vertex projective 
ambiguity, and y and x' belong to a same site, so y^^ • x' does not contain the site ambiguity. 
Finally we have to provide the projective invariance of Jw with respect to the start and 
end points of each simply connected subwalk. In our case of the local Weyl algebrae this 
invariance is obvious, because of elements x • y— 1 for different vertices commute. 

With the structure of the walks introduced, the simple analysis of the determinant gives 
immediately 

J.,/3 = E (-)'^(^-^) • Jw„,' (3-54) 
all 

where the sum is taken over all the walks of the homotopy class a A + (3 B given and the 
system of the outlets of the walks fixed. 



3.6 Monodromy operator 

Consider now another interpretation of the two dimensional kagome lattice. 

Let now to each vertex of the lattice the local L-operator is assigned. Instead of lo ( 2.931 ) 
in the definition of L, use 

x' = C"^ • X , y' = C • y , (3.55) 

where 

C = C{x-^u,y-^\N) (3.56) 

and 

C(u,w) = - u^^w + KW, (3.57) 

For the A-type triangle P let the in - edge variables be xp, yp, zp so that out edge 
variables are Xqp, y^p and Zcp, c = a^^ b. These notations are shown in Fig. Surely 
LYBE for L-operators means that for xp, yp, zp given, x^p, yf,p and Zcp are the same for 
the right hand side YBE graph V- 

Consider now the whole toroidal kagome lattice. We are going to assign x, y and z to 
some minimal set of the edges so that the variables of all other edges can be restored via 
functions Cj^p. 

To do this, cut the torus along some line, shown as the dashed line in Fig. ^ Call 
this line 'the string'. The edge variables along the string we'll denote as Xj, y^ and Zj. It is 
useful to draw the string so that it intersects all x and z lines once, and y-lines twice (i.e. 
the homotopy class of the string is ±(2 ^ — ;B), an orientation of the string and so a sign are 
unessential) Note, Xj, y^ and Zj introduced we assign to the edges which are right-touched 
to the string. 
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Now switch on the L-operator game with the edge variables. We interpret it as the shift 
of the string. We enumerate the Hnes so that the triangle P = Pa is surrounded by 

the lines , and z^+j , so as the L - operators are 

Lx.,y^.({u,w}i^a,t.pJ , Lx.,zJ{u,w}2,a'=-h.Po) , Ly^,zJ{u, w}3_<j, bfc-. pj . (3.58) 

The L-operator game allows us to restore all the edge variables for the lattice, including the 
left-touched to the original string variables x^ , y^ , z^ . 

Thus for given variables from the right side of the string we obtain the analogous val- 
ues from the left side of the string as functionals of the given variables. Thus the map 
corresponding to the kagome lattice and the string chosen appears: 

T(L) : {x„y,,z,} ^ {x:,y',z',}. (3.59) 



As the operator, T(L) is ordered product of all L ( 3.58 ). Define A < B if the ordered 
product of A and B is A ■ B. Then in r(L) 

Ly^.,z,+, < Lx,,z.+, < Lx,,y^ , (3.60) 
Lx,,y^ < Lx,,z,+,+i < Ly^..z,+,+i • (3.61) 

These relations are enough to restore T(L). Operator r(L) resembles the monodromy matrix 
in 2D. The difference is that instead of the distinguished point in 2D monodromy matrix 
(i.e. the point where the transfer matrix is torn), in 31? we have the distinguished string. 
Now, what should stand for a "trace" of T(L). Consider the system 

Xj = Xj , y'j = , z^- = Zj (3.62) 

on some left module element (pQ. Here are 3 M equations, 3 Af — 1 from them are independent 
due to 

n X, z, = n 4 y; 4 ' (3-63) 

3 3 

where it is implied that all xp,yp,zp VP arc commutative, this is the c onse quence of the 
commutability of primary Xj , y^ , Zj . Then solve 3 M — 2 equations of ( [3.62 ) leaving two 
variables, up to unessential signs and powers of q : 

A = l[ qy^z^, B = n ^'^3- (3-64) 

j j 

A single equation rests for A and B, and amusingly this equation coincides with the quantum 
determinant relation (j)^ ■ i{A,B) = 0. So in this sense ^{A,B) = is the trace of the 
monodromy operator. 

Note however, ^{A,B) is the invariant curve, this was established in the previous 
section, so it is not necessary to consider (pQ ■ ^{A,B) = 0. The actual problem for the 
further investigations is to diagonalize J ( A, B) for A and B arbitrary. 

4 Discussion 

Conclude this paper by an overview of the problems to be solved and the aims to be reached. 
The approach proposed gives a way for their solution. 
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First, mention the problems of the classification of the map 



R : {a„b„c„d,} ^ {a;,b;.,c;,d;.}, J = 1,2,3, (4.1) 

in general. The aim is to classify all conserving symplectic structures of the body B. We 
have discussed only the local case, when the variables, assigned to different vertices com- 
mute, and the scalars (spectral parameters) are conserved. We suppose, such case is not 
unique, and there are another ways to remove the projective ambiguity. The simplest case 
to be investigated is to consider all the variables a, b, c, d for each vertex as matrices with, 
for example, non-commutative entries, but with this entries commutative for any two ver- 
tices. The matrix structure may be common for all vertices, and thus we would have no 
commutation between different vertices in general. Another simple possibility is another 
kind of locality, the case when the dynamical variables commute while do not belong to a 
same site. 

Note, once our locality is imposed, the Weyl structure appears immediately. Thus the 
Weyl algebra is the consequence of the locality technically, but a principal origin of the Weyl 
algebra is mysterious. 

Next fundamental problem is the quantization of Korepanov's matrix model mentioned 



above. The conservation of complete algebra X, ( 2.76 , 2.77 ), means that we can not use 
( 2.79 ) to fix all the ambiguity of Korepanov's equation. Analysis of ( ^.74 ) plus some extra 



(but natural) symmetry conditions allows to fix the functional map r : Xj i-^ Xj up to one 
unknown function of three variables. The problem of the Tetrahedron equation for these r 
is open. All these are a subject of a separate paper. 

Pure technical problem to be mentioned is the investigation of q - hypergeometrical 



function a, eqs. (2.53,2.5?:) 



Another interesting thing is functional L - operators and LYBE related to them. The 
map given by L, eqs. (2.93, ^.94 ), is a bi-rational one. Note, the case of linear L coincides 



with Korepanov's X. Thus the rational case of it as well as the general case of the bi-rational 
transformation have good perspectives for the investigation. 

The main set of problems for further investigations is connected with the integrals of U . 
^{A, B) seems to be not constructive. The aim is to calculate the spectrum of it, and to 
calculate U as a function of its integrals. Possible approach is functional equations for the 
integrals of motion, that should follow from the determinant or topological representation 
of J. 

Another possibility is a way resembling the Bethe ansatz in 2D should exist in 3D, i.e. 
a way of a triangulation of U with a help of some artificial operators. If such way exists, it 
must based on the linear problem derived. 
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